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Abstract 

This paper describe a numerical simulation method for the interaction between 
laser pulses and low density plasmas based on hydrodynamic approximation. We 
investigate Backward Raman Amplifying (BRA) experiments and their variants. 
The numerical results are in good agreement with experiments. 
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1 introduction 



Recently, the concept of laser amplifying by backward raman scattering(Backward 
Raman Amplifiers, BRA) in plasma was presented [1],[2],[3],[4]. In which a 
short seed pulse is amplified by a counter-propagating long pumping pulse 
when the resonance relation Ugeed + ^plasma = '^pump is satisfied. In the non- 
linear regime, the seed pulse is strongly amplified and compressed temporally. 
Theoretically, unfocused intensities of IQ^'^W / cm? in 50 fs pulses are accessible 
by this technique. 

Up to now two schemes of BRA have been presented: the amplifier in non- 
linear regime and transient regime. The basic concepts of the first scheme are 
based on the vr-pulse solutions of laser-plasma coupled equations [2]. When the 
amplitudes of the plasma wave and the seed pulse are both quite large, the 
pumping pulse's energy is absorbed by the front of the seed pulse. So the front 
edge of the seed pulse is amplified while the back edge will be reduced for the 
energy feedback into the exhausted pumping pulse. In this way, the seed pulse 
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is amplified and compressed. Finally, the seed pulse evolves to a narrow wave 
train which is called a tt pulse. 

The scheme of transient BRA is based on the matching of the plasma's density 
gradient and the pumping pulse's frequency chirp [4] [5]. In the amplifying 
process, the strict matching of the two pulses and the plasma frequency is 
always needed. If we let the pumping pulse chirp and the plasma density have 
a gradient, and set the plasma frequency gradient just be a half of pumping 
chirp gradient, the seed pulse will be amplified only at a fixed position of 
itself. By passing a plasma of sufficient length, the amplifying take place and 
a ultra-fast laser pulse formed. 

Each scheme is very difficult to be realized. For the limit of plasma technology, 
the density of gas-jet plasma is limited for about lO^^cm"^ now. To get distinct 
amplification, the seed pulse must travel several millimeters in plasma. The 
amplification require rigorous match of the frequencies. Keeping the uniformity 
of the plasma is a great challenge in experiments. In the transient scheme, we 
need even longer plasma length to compensate the lower efficiency. 

However, some good results of amplifying experiments have been posted recently[6] [7] . 
Especially, the experiments informed that the non-linear BRA scheme was pos- 
sible. In the experiments by W.Cheng et al, the seed pulses were amplified to 
one order higher than pumping pulse and the amplifying ratio was several 
thousands. 

Consequently, the numerical simulation of BRA becomes a valuable problem in 
computational plasma physics. However, compared with the common plasma 
simulations there are some special difficulties because of the long-distance 
travelling wave[8][9][10]. First,thc computing time and memory needed is very 
large. Second, the investigation requires low-error simulation technology for 
the low increase rates of Raman effects . For example, we need a great number 
of particles to get rid of the numerical wave-breaking and numerical heating in 
Particle-In-Cell method[ll]. Besides, we need a high-order difference scheme 
of the electromagnetic ffeld equation to overcome the difficulties of numeri- 
cal dispersion while applying it to the long-distance traveling wave. So the 
simulations of BRA are mainly focused on two regime now: 1) averaging the 
effects of electromagnetic fields to reduce the space-time steps, such as the 
APIC method[12]. 2)using the amplitudes coupled equations [13], [14], [15] or 
linearized plasma wave equations to get approximately analysis. 

An attracting idea is applying a moving window which following the seed pulse 
to the computation. This technology is convient to investigate the interactions 
between laser and low-density plasma. In the corona region, the plasma density 
is about Yq '^o The group velocity of laser pulse is very close to light speed. 
So the calculated window is shorter than ^ of the total interaction region. 
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There is a simulation method which apphes a moving window on the stan- 
dard PIC simulation (Moving window PIC Simulation, MWPIC) [16]. However 
MWPIC needs the space-time step match((ix = cdt) to avoid the interpolation 
error. The step setting is unstable at 2-dimensions and 3-dimensions explicit 
differential schemes. That means we need some implicit differential schemes. 
The latter technology are still developing now. In addition, MWPIC need a 
very large number of particles to get acceptable precision too. 

For the above reasons, we developed a numerical simulation technique based on 
hydrodynamic approximation. In this technique, we only use the cold-plasma 
approximation which is often used in theoretical researches. The numerical 
trials indicate that the technique is suitable for the traveling wave problems 
in plasmas and can give results in concordance with experiments. Besides, the 
technique has a clear physical picture and the results can be easily analyzed. 

We will discuss the BRA equations and the simulation technology in section 
2 and section 3, then analyse the numerical results in section 4. 



2 the basic equation of BRA 



In the typical instance of BRA, we can introduce the cold-plasma approxima- 
tion and the transverse canonical momentum conservation [171: v = u + = 

L J 7mc 

— * 

u + ca, where u is the longitude velocity and A is the vector potential at 
Coulomb gauge while a is the normalized vector potential: a — Because 
the electrons' velocities are far less than the light speed,we can ignore the 
relativistic factor 7. 

The plasma hydrodynamic equation can be written as: (We ignore the thermal 
pressure here. We can recover it to the rhs. of the second equation whenever 
needed) 



dn 
'di 

dv 

di 



eE e{v X B) 



m 



mc 



(1) 
(2) 



Supposing the plasma flow are non-vertex before the laser action and applying 
the cold-plasma approximation,we get the plasma hydrodynamic equation: 
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+ (/3 . V)n = -n( V • /3) - (a ■ V)n (3) 



d{ct) 



1^ + (/3 • V)/3= V0 - ^Va^ - V(/3 • a) (4) 
The electromagnetic field satisfies the wave equation and Poisson equation: 



'n - Hi) (6) 



Applying the coordinate transformation 

C^x-ct (7) 
V = t (8) 

we get the wave equation in the new variables: 



)^--^-r^ (9) 



and the Poisson equation has no change. 

To solve the wave equation numerically,we introduce the medium variable t/^: 

c d-q 

Finally,we get a pair of first order differential equations: 



9V' Anne^ - , V72 - 
^ = + V| a 



With the same transformation,the hydrodynamic equations became 
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(11) 



The equations (6), (10), (11) are the basic equations to describe the laser- 
plasma interactions in BRA. 



3 numerical simulation method 



The electromagnetic equations (6), (10) describe the electromagnetic wave 
transmitting in plasma. In BRA scheme, the transverse dynamic effects arc 
mainly slow processes such as self-focusing. We can use the periodic boundary 
conditions in the transverse directions, apply FFT calculation to eliminate the 
transverse derivatives and use direct elimination technology in the longitudinal 
direction. So we will focus on the 1-dimension scheme at following discussion. 
However, the generalization to multi-dimensions is straightforward. 

The elimination of Poisson equation is quite simple. In our moving window, 
the right boundary of calculating region are always at the front of seed pulse. 
We ignore the interactions between the pumping pulse and the origin plasmas 
for the uniform of the pumping pulse, so all laser-plasma interaction can not 
affect the region in the right side of the boundary. Then we can assume that 
there are no electrostatic fields in the right region of seed pulse: ^ = then 
— = ... = const. The differential scheme is 

47r6^ 

4-1 = ^5x^{ni - Ni) - + 24 (12) 



To get the numerical solution of the wave equation, we use a generation of 
R.Liu's leapfrog-upwind scheme[18]: 



5t 



-2c 



5x 
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Where 



K= 2c|^ (13) 
dx 

m= ^(5 + 3k -2^2) (14) 

«r'=<+i'-(l-2m«:) «+!-<) 

+(1 - m)fi;«+2 - <+i) + (1 - m)fi;« - <_i) 

+c5t(C+i+0 (15) 



The hydrodynamic equation (11) can be solved with leap-frog method too. 
However,we must apply the FCT technology [19] to avoid the numerical oscil- 
lation while introduce little numerical damping. The Flux Corrects are used 
in all leap-frog steps by rewritten to convection equations: 



F means the convect terms and S means source terms,then the leap-frog 
scheme is: 



ur' = ur' - 2^(^+1/2 - ^-1/2) + 2At^r (17) 



To use FCT technology, first we calculate the low-order flux of plasma den- 
sity(there is no source term in the continuum equation): 



^t') (18) 

(19) 

(20) 
(21) 
(22) 

7/ is a coefficient about ^ to |. The smaller rj means less error but the more 
numerical oscillation. Then the high-order flux is: 



fS = ^(/^+i + - ^(^2 + /.-i) (23) 
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We don't introduce the correct of the velocity's flux because it wiU cause much 
damping in the simulation. So the velocity flux can use the same scheme as 
the high-order density flux. 

However, the difference scheme of plasma density does not preserve positivity. 
To get rid of the difficulties of negative electron density, we apply a hmiter of 
density flux: the out-flux density can not exceed the present electron density 
of each point. The physics means of this limiter are clear. 



4 numerical results and the analysis 

W.Cheng's experiment results [7] were the best experiment validate of BRA 
scheme. In this experiment, the author amplified the seed pulse in a gas- 
jet plasma. In the plasma with 2mm effective length, the intensities of seed 
pulses were amplified several thousands times and the length of seed pulses was 
compressed from 550fs (origin) to about 150fs. The recorded intensity reached 
1.7 X lO^^W/cm^ which was higher than pump intensity by more than 1 order 
in magnitude. The author recorded the broadened ( 800 fs ) and compressed 
(< 200/s) output pulses both in the experiment. Which informed that the 
amplifying process had been in non-linear region. 

We investigated the process numerically. As the experiments setting up, we set 
the seed pulse length to 550 fs and set the pumping pulse length to 10 ps. The 
plasma's length was about 2 mm and the plasma's density was 1.1 * 10^^ cm~^ 
. The wavelength of pumping pulse was 0.8 /im and the wavelength of seed 
pulse was set to matching with the plasma and pumping pulse. 

At the peak top of the pumping pulse, we let the two pulses collide each other. 
For the uncertainty of the collide points and the plasma temperatures, the 
simulating results varied a httle. When the interaction finished, the seed pulse 
evolved to narrow peek which FWHM about 150/s — 170/s and peak intensity 
about 1.5 * lO^^W/cm"^ — 1.9 * lO^^W/cm'^ . The typical simulate results were 
shown at figure 1. We can find that the vr-pulse shape easily. Which means 
the non-linear amplifying region was reached. Considered the uncertainly of 
the plasma temperature and length, we can find good agreement between the 
simulations and the experiments. 

From the simulation results, we can find the evolving process of seed pulse. Let 
the seed pulse collide with different length pumping pulse,we can image the 
evolving of the seed pulse: it was broadened to about 750 fs and compressed 
later. At the end, it evolved to a pulse shorter than 160 fs. 

The typical results are shown at figure 2. 
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figure l:The output intensity versus ^ after amplified. 




figure 2: The amplified pulse shape of the seed by different length pumping 

pulse. 

The simulating process shows the non-linear character of plasma wave. When 
the amplitude of plasma wave is quite large, the wave's shape are different from 
sine wave, the typical wave form are shown in figure 3. As we see, the difference 
are quite evidence in our results. In fact, the effective 5n/n can exceed 1 when 
the amplifier reaches the exhaust region. This phenomena made an adding 
stability of amplifier from wave-breaking. 

However, our results show that when the amplifying length reach to 2mm, the 
plasma wave amplitude was close to cold plasma wave-breaking amplitude. 
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figure 3: The typical wave-shape of density oscillating when the amplitudes 

become large. 

Which means, if wc want to make an practical amplifier ,we must increase 
the plasma density to get higher stabilities. Although It was very difficult to 
get dense and uniform plasma in experiments, we simulated the BRA process 
in a quite dense plasma {n — O.OAricr)- The results were coarse in precision 
but show that the amplifying processes were much more stable and easier to 
reach the nonlinear-regime in a more dense plasma. The higher increase rates 
of raman effects and higher stability compensated the lower conversion effects 
in the more dense plasmas. The results of our simulation are shown in figure 
4. In this density, we need shorter than 1mm plasmas to generate obviously 
amplified intensities. It means that the increasing of uniform plasma density 
should be the most important improvement at future experiments. 

Our method can be easily improved to include the effects as pumping heat, plasma 
un-uniformity, inverse-Bremsstrahlung, or even linear landau damping, tunnel- 
ionized (by phenomena) effects et al. The investigation of these effects in BRA 
are proceeding now. 



5 conclusion 

We introduced a hydrodynamic simulation method which was suitable for the 

interaction of lasers and corona plasmas. The method is convenient for the 
simulation of ultra-fast laser's transmitting. We simulated the Backward Ra- 
man Amplifiers experiments numerically, and the results show good agreement 



2568 



0.050 -I 



to 

3 0.045 - 




5000 10000 15000 20000 25000 

5=x-ct by plasma units 

figure 4: The pulse shape of the seed after being amplified by different length 
pumping pulse;the unit is the natural of plasma. 

with experiments. The following study will appear in another paper. 
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